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Abstract 



We establish a connection between the strong solution to the spa- 

lO ■ tially periodic Navier-Stokes equations and a solution to a system 

^ , of forward-backward stochastic differential equations (FBSDEs) on 

^vq I the group of volume-preserving diffeomorphisms of a flat torus. We 

^J ■ construct representations of the strong solution to the Navier-Stokes 

equations in terms of diffusion processes. 

o 

O ■ Keywords: Navier-Stokes equations, forward-backward SDEs, diffeomor- 

phims group 

c^ ■ 1. Introduction 

The classical Navier-Stokes equations read as follows: 

d 

— M(t, x) = — (m, V)u{t, x) + uAuit, x) — Vp(t, x), 

divu = 0, ^^^ 



u{0,x) = —uo{x] 



where uo{x) is a divergence- free smooth vector field. We fix a time interval 
[0,T], and rewrite equations (1) with respect to the function 

u{t, x) = —u{T — t,x). 

Problem (1) is equivalent to the following: 

d 

Tru{t, x) = — (n, V)n(t, x) — z/An(t, x) — Vp{t, x), 

divM = 0, (^) 

u{T,x) = uo{x), 

where p{t, x) = p{T — t,x). 

In what follows, system (2) will be referred to as the backward Navier- 
Stokes equations. To this system we associate a certain system of forward- 
backward stochastic differential equations on the group of volume-preserving 
diffeomorphisms of a fiat torus. For simplicity, we work in two dimensions. 
However, the generalization of most of the results to the case of n dimensions 
is straightforward. The necessary constructions and non-straightforward gen- 
eralizations related to the n- dimensional case are considered in the appendix. 

Assuming the existence of a solution of (2) with the final data in the 
Sobolev space H°' for sufficiently large a, we construct a solution of the 
associated system of FBSDEs. Conversely, if we assume that a solution of the 
system of FBSDEs exists, then the solution of the Navier-Stokes equations 
can be obtained from the solution of the FBSDEs. In fact, the constructed 
FBSDEs on the group of volume-preserving diffeomorphisms can be regarded 
as an alternative object to the Navier-Stokes equations for studying the 
properties of the latter. 

The connection between forward-backward SDEs and quasi-linear PDEs 
in finite dimensions has been studied by many authors, for example in [9], 
[18], and [21]. 

Our construction uses the approach originating in the work of Arnold [1] 
which states that the motion of a perfect fiuid can be described in terms of 
geodesies on the group of volume-preserving diffeomorphisms of a compact 
manifold. The necessary differential-geometric structures were developed in 
later work by Ebin and Marsden [10]. We note here that [1] and [10] deal 
only with differential geometry on the group of maps without involving prob- 
ability. 



The associated system of FBSDEs is solved using the existence of a so- 
lution to (2), and by applying results from the works of Gliklikh ([13], [14], 
[15], [16]). The latter works use, in turn, the approach to stochastic differen- 
tial equations on Banach manifolds developed by Dalecky and Belopolskaya 
[4], and started by McKean [19]. Conversely, a solution of (2) is obtained 
using the existence of a solution to the associated FBSDEs as well as some 
ideas and constructions from [9]. However, unlike [9], we work in an infinite- 
dimensional setting. 

Representations of the Navier-Stokes velocity field as a drift of a diffu- 
sion process were initiated in [23] and [24]. A different system of stochastic 
equations (but not a system of two SDKs) associated to the Navier-Stokes 
system was introduced and studied in [5]. This system also includes an SDE 
on the group of volume-preserving diffeomorphisms, but is not a system of 
forward-backward SDKs. Also, we mention here the works [2] and [3] dis- 
cussing probabilistic representations of solutions to the Navier-Stokes equa- 
tions, and the work [6] establishing a stochastic variational principle for the 
Navier-Stokes equations. Different probabilistic representations of the solu- 
tion to the Navier-Stokes equations were studied for example in [17] and [7]. 
We note that the list of literature on probabilistic approaches to the Navier- 
Stokes equations as well as connections between finite-dimensional FBSDEs 
and PDEs cited in this paper is by no means complete. 

The method of applying infinite-dimensional forward-backward SDEs 
in connection to the Navier-Stokes equations is employed, to the authors' 
knowledge, for the first time. 

2. Geometry of the difFeomorphism group of the 2D 
torus 

Let T^ = S^ X S^ be the two-dimensional torus, and let if°(T^), a > 2, 
be the space of if"-Sobolev maps T^ — )• T^. By G" we denote the subset of 
iJ"(T^) whose elements are C^-diffeomorphisms. Let G^ be the subgroup of 
G" consisting of diffeomorphisms preserving the volume measure on T^. 

Lemma 1. Let g be an H°'-map and a local dijjeoniorphism of a finite- 
dimensional compact manifold M , F be an H"-section of the tangent bundle 
TM . Then, F o g is an H'^-map. 

Proof. See [15] (p. 139) or [10] (p. 108). D 



Let Rg denote the right translation on G°, i.e. Rg{r]) = r] o g. 

Lemma 2. The map Rg is C°° -smooth for every g G G". Furthermore, for 
every rj G G", the tangent map TRg restricted to the tangent space Tj^C^ is 
defined by the formula: 

TRg : TrjG" -> T^ogG", X^Xog. 

Proof. The proof easily follows from the a-lemma (see [10], [15], [16]). D 

Lemma 3. The groups 0°" and G" are infinite- dimensional Hilbert manifolds. 
The group Gy is a subgroup and a smooth submanifold ofG°'. 

Lemma 4. The tangent space TgG°' is formed by all H"' -vector fields on T^. 
The tangent space T^G^ is formed by all divergence-free H"" -vector fields on 

The proof of Lemmas 3 and 4 can be found for example in [10], [15], [16]. 

Lemma 5. Let X G TgG" be an H'^ -vector field on T^. Then the vector field 
X on G°' defined by X{g) = X o g is right-invariant. Furthermore, X is 
C'^ -smooth if and only if X & H^^^ . 

Proof. The first statement follows from Lemma 2. The proof of the second 
statement can be found in [10]. D 

The vector field X on G° defined in Lemma 5 will be referred to below 
as the right-invariant vector field generated by X G T^.G'^. 

Let g G G", X,Y & TgG"'. Consider the weak (-, ■ )o and the strong 
( ■ , ■)a Riemannian metrics on G" (see [16]): 

{X{g), Y{g))o = f {X o g{e), Y o g{e))de, (3) 

{X{g),Y{g))^= [ iXog{9),Yog{e))de 

JT2 



+ / i{d + 5)''Xog{e),{d + 5rYog{e))de (4) 

JT2 

where d is the differential, 6 is the codifferential, X and Y are the right- 
invariant vector fields on G° generated by the if "-vector fields X and Y. 
Metric (3) gives rise to the L2-topology on the tangent spaces of G", and 



metric (4) gives rise to the if "-topology on the tangent spaces of G" (see 
[16]). li g E G^, then scalar products (3) and (4) do not depend on g. 
Moreover, for the strong metric on G", we have the following formula: 

{X{g),Y{g))^= f {Xog{e),{l + ArYog{e))de 

where A = {d6 + 6d) is the Laplace-de Rham operator (see [22]). 
Let us introduce the notation: 

Z+ = {{ki, fcs) eZ^ ■.ki>OoT ki = 0,k2> 0}; 



k = {ki, k^) e Z+, k= {k2, -ki), \k\ = ^Jkl + kl k-9 = ki9^ + M2, 



and the vectors 



1 ., .. / A;2 A r. rn-. 1 • n m / h 



MO) = J^^^-os{k . 6) ^_IJ , B,{e) = j^;^,Mk ■ 9) ^_^^ 

Let {^fc(5'), -BA;(5')}fcgz+u|o} ^^ ^^^ right-invariant vector fields on G°' gener- 
ated by {Afc,5fc}fcg^+u^o|, i.e. 



a 



Ak{g) = Akog^ Bk{g) = Bkog, geC 

Aq = Aq, Bq = Bq. 

By w- lemma (see [15]), A^ and B^ are C°°-smooth vector fields on G". 

Lemma 6. The vectors Ak{g), B^ig), k e Z^ U {0}, g e G", form an 
orthogonal basis of the tangent space TgG^ with respect to both the weak 
and the strong inner products in TgG^ . In particular, the vectors A^, B^., 
k G Z^U{0}, form an orthogonal basis of the tangent space T^G^. Moreover, 
the weak and the strong norms of the basis vectors are bounded by the same 
constant. 

Proof. It suffices to prove the lemma for the strong norm. Let us compute 
A^Ak. Note that the vectors jg and t|t form an orthonormal basis of M?. 



Let us observe that by the identity {k, V) cos(/c ■ 9) = 0, 6Ak = 0. Hence 
d6 Ak = which imphes AAk = SdAk. We obtain: 

- 1 k 

^'"|iF "''*'"* 1*1' 

- - 1 k - 

AAk = 5dAk = Trn^cos{k-e) — = \k\^Ak, 

\ \ \ \ 

A" Ak = |A;r cos{k ■ 0) -^ = |A;|^"Afc. 

This and the volume-preserving property of (7 G G" imply that 

(B^ig), Akig))a = (^^, ^)a = (1 + |A;p")(5^, ^)l, = 0, 
11^.(5)11^ = ll^ll^ = (1 + \kn P,||L= 2^' d^r'" + 1) 

where || ■ ||q, is the norm corresponding to the scalar product ( ■ , ■)a- Thus, 
27r2 ^ ||Afc(5')||^ ^ 47r2. Clearly, for the ||Sfc(5f)||^ we obtain the same. D 

It has been shown, for example, in [10] and [16] that the weak Riemannian 
metric has the Levi-Civita connection, geodesies, the exponential map, and 
the spray. Let V and V denote the covariant derivatives of the Levi-Civita 
connection of the weak Riemannian metric (3) on G° and G^, respectively. 
In [10] (see also [16], [15]), it has been shown that 

V = P o V 

where P : TG" — ?■ TG" is defined in the following way: on each tangent 
space TgG", P = Pg where Pg = TRg o P^ o TRg^i, TRg and TRg-i are 
tangent maps, and Pg : TgG" — )■ TgGy is the projector defined by the Hodge 
decomposition. 

Lemma 7. Let U he the right-invariant vector field on G°' generated by an 
H""^^ -vector field U on T^, and let V he the right-invariant vector field on G°' 
generated hy an H"' -vector field V on T^. Then Vyll is the right-invariant 
vector field on G" generated hy the H°' -vector field VyU on T^. 



Lemma 8. Let U he the right-invariant vector field on Gy generated by a 
divergence-free H"^^^ -vector field U on T^, and let V he the right-invariant 
vector field on G° generated hy a divergence- free H"" -vector field V on 
T^. Then Vytl is the right-invariant vector field on G^ generated hy the 
divergence-free H'^ -vector field PeVyU on T^. 

The proofs of Lemmas 7 and 8 follow from the right-invariance of covariant 
derivatives on G° and Gy (see [16]). 

Remark 1. The basis {A^, -Bfc}fcgg+u|o} of T^G^ can be extended to a basis 
of the entire tangent space T^G"". Indeed, let us introduce the vectors: 

MO) = ^ cos{k . 9) (^j) , Bm = ^ sin(fc . 9) (^j) , fc G Z^ 

The system A^, Bk, k G Z^ U {0}, Ak, Bk, k G Zj', form an orthogonal basis 
of TgG". Further let Ak and Bk denote the right-invariant vector fields on G" 
generated by Ak and Bk- 

3. The FBSDEs on the group of difFeomorphisms of the 
2D torus 

Let /i : T^ — )■ M^ be a divergence- free iJ°"''^-vector field on T^, and let h be 
the right-invariant vector field on G" generated by h. Further let the function 
V{s, ■ ) be such that there exists a function p : [t, T] — )■ H°'^^(T'^, M) satisfying 
V{s, ■ ) = Vp(s, ■ ) for all s G [t,T]. For each s G [t,T], V{s, ■ ) denotes the 
right-invariant vector field on G" generated by V{s, ■ ) G if°'(T^,]R^). 

Let E' be a Euclidean space spanned on an orthonormal, relative to the 
scalar product in E, system of vectors {e^, ef , e^, e^}fcgz+ \k\^N- Consider the 
map 

a{g)= Yl Ak{g)®et + Bk{g)®e^.geG\ 

fcez^ujo}, 

i.e. cr((y') is a linear operator E — )■ T^G" for each g G G". 

Let (f2, J^, P) be a probability space, and Ws, s G [t, T], be an i?-valued 
Brownian motion: 



Ws= Y. il^kis)et + Pi{s)e^) 



fcez+u{o}, 
\k\<:N 



where {/3^, l3k}kez+u{o},\k\<:N i^ a sequence of independent Brownian motions. 
We consider the following system of forward and backward SDEs: 

dZ'/ = Y^^'^ds + ea{Z'/)dWs, 

dY'^'^ = -V{s, Zl'^)ds + Xl'^dWs, (5) 

Z^ = e; Yj, = n[Zr^ ). 

The forward SDE of (5) is an SDE on G" where G^ is considered as a 
Hilbert manifold. Stochastic differentials and stochastic differential equations 
on Hilbert manifolds are understood in the sense of Dalecky and Belopol- 
skaya's approach (see [4]). More precisely, we use the results from [15] which 
interprets the latter approach for the particular case of SDEs on Hilbert man- 
ifolds. The stochastic integral in the forward SDE can be explicitly written 
as follows: 

r a{Z'/)dWr = Y. f MKnd/3tir) + B^iZ^nd/^iir). (6) 

Let us consider the backward SDE: 

Y^''' = h{Z'/) + V{r,Z'/)dr- X'/dW^. (7) 

J s J s 

Note that the processes V{s, Z^'^) = V{s, ■ ) o Z*'*^ and h{Zr^^) = h o Zrf are 
if "-maps by Lemma 1. Therefore, it makes sense to understand SDE (7) as an 
SDE in the Hilbert space //"(T^R^). Let J=, = a{Wr,r e [0,s]). We would 
like to find an J^^-adapted triple of stochastic processes (2'*'^, Y^'"^, X^''^) solv- 
ing FBSDEs (5) in the following sense: at each time s, the process (2'*'*^, Y^''^) 
takes values in an if "-section of the tangent bundle TG^. Namely, for each 
s e [t,T] and w G fi, Z*'" G G°, F/'*^ G T^tM^. Therefore, the forward 
SDE is well-posed on both G" and G" , and can be written in the Dalecky- 
Belopolskaya form: 

dZ'f = exp^t.ejr/'^ds + ea{Zl'')dWs} 

or 
dZ'/ = exp^t.ejF/'^ds + ea{Zl'^)dWs} 

where exp and exp are the exponential maps of the Levi-Civita connection 
of the weak Riemannian metrics (3) on G°' and resp. G" . Below, we will 



show that using either of these representations leads to the same solution of 
FBSDEs (5). 

Finally, the process X*'*^ takes values in the space of linear operators 

/:(E,i7°(T2,M2)),i.e. 

Xl'^= Y. X's^^^ + X's^'^e^, (8) 

fcGZ+u{0},|fc|s:Af 

where the processes X^^ and X^^ take values in if"(T^,M^). 

Remark 2. The results obtained below also work in the situation when 
the Brownian motion Ws is infinite dimensional (as in [8]). Namely, when 
Ws = Efcez+ujo} ^k(3k ® e^ + hf^k ® ef where a^, bk, fc G Z+ U {0}, are 
real numbers satisfying J2kGZ+um I'^fcl^ + 1^^=!^ ^ '^^ However, this requires 
an additional analysis on the solvability of the forward SDE based on the 
approach of Dalecky and Belopolskaya [4] since the results of Gliklikh ([13], 
[15], [16]) applied below are obtained for the case of a finite-dimensional 
Brownian motion. 



4. Constructing a solution of the FBSDEs 

4.1 The forward SDE 

Let us consider the backward Navier-Stokes equations in M^: 

y(s, 9) = h{9) + ^ [Vp(r, 9) + {y{r, 9), V)y{r, 9) + uAy{r, 9)] dr, 
dwy{s,9) = 
where s G [t,T], 9 ^ T^, A and V are the Laplacian and the gradient. 

Assumption 1. Let us assume that on the interval [t,T] there exists a solu- 
tion {y{s, ■ ),p{s, ■ )) to (9) such that the functions p : [t,T] -> H°'+'^(T'^,R) 
and y : [t, T] — )■ H°'~^^(T'^,M.'^) are continuous. 

Clearly, y{s, ■ ) G T^G^. Let {F,*-^^, ^"/'''■^jfcez+ulol ^e the coordinates of 



ru{0} 
y{s, ■) with respect to the basis {Afc, -B/cj^g^+^jro}, i-e. 

yis,9)= Yl ys''''Md) + Y'*''M9). 
fcezju{0} 



Let Ys{-) denote the right-invariant vector field on G" generated by the 
solution y{s, ■), i.e. Ys{g) = y{s, ■) o g. On each tangent space TgG"', the 
vector Ys{g) can be represented by a series converging in the //"-topology: 

ys{g)= Yl Ys''^M9) + Y:-^'''B,ig). (10) 

fcez+u{o} 

In this paragraph we will study the SDE: 

dZ'/ = Ys{Zl'')ds + ea{ZY)dWs. (11) 

Later, in Theorem 6, we will show that the solution Z^ to (11) and the 
process y/'*^ = Ys{Zl'^) are the first two processes in the triple (^5'^, F^*'^, -^i'*^) 
that solves FBSDEs (5). 

Theorem 1. There exists a unique strong solution Z*'"^, s G \t,T], to (11) 
on G^ , with the initial condition Z^^ = e. 

Proof. Below, we verify the assumptions of Theorem 13.5 of [16]. The latter 
theorem will imply the existence and uniqueness of the strong solution to 
(11). Note that, if sum (6) representing the stochastic integral /^ a{Zl''^) dWg 
contains only the terms Aq{I3q{s) — Poif)) and Bq{I3q{s) —(3^(1:)), i.e., infor- 
mally speaking, if the Brownian motion runs only along the constant vectors 
Ao and Bq, then the statement of the theorem follows from Theorem 28.3 of 
[16]. If sum (6) contains also terms with A^ and B^, k ^ Z^, or, informally, 
when the Brownian motion runs also along non-constant vectors Ak and Bk, 
k e Z2, then the assumptions of Theorem 13.5 of [16] require the bounded- 
ness of Ak and Bk with respect to the strong norm. The latter fact holds by 
Lemma 6. 

Hence, all the assumptions of Theorem 13.5 of [16] are satisfied. In- 
deed, the proof of Theorem 28.3 of [16] shows that the Levi-Civita con- 
nection of the weak Riemannian metric (3) on G" is compatible (see Def- 
inition 13.7 of [16]) with the strong Riemannian metric (4). The function 
(^(9) = J2 kez+ u{o}, \k\^N ^^kig) ® e^ + Bk{g) ef is C~-smooth since Ak and 
Bk are C°°-smooth. Moreover, by Lemma 6, a{g) is bounded on G^. Next, 
since y : [t,T] — )■ //""^""^(T^, M^) is continuous, then it is also bounded with 
respect to (at least) the i/"-norm. Hence, the generated right-invariant vec- 
tor field Ys{g) is bounded in s with respect to the strong metric (4), and 
it is at least C^-smooth in g. The boundedness of Yg in g follows from the 
volume-preserving property of g. D 

10 



Theorem 2. There exists a unique strong solution Z\'^ , s G [t,T], to (11) 
on G"', with the initial com 
solution to SDE ill) on G^ 



on G"', with the initial condition Z^'^ = e. This solution coincides with the 



Proof. Consider the identical embedding i : G^ — )■ G". By results of [4] 
(Proposition 1.3, p. 146; see also [16], p. 64), the stochastic process i{Zl'^) = 
Zl''^, s G [t,T], is a solution to SDE (11) on G", i.e. with respect to the 
exponential map exp. This easily follows from the fact that Ti : TG^ — )■ TG", 
where T is the tangent map, is the identical imdedding, and that z(exp(X)) = 
exp(TzoX). The solution Z*'^ to (11) on G" is unique. This follows from the 
uniqueness theorem for SDE (11) considered on the manifold G" equipped 
with the weak Riemannian metric. Indeed, a{g) and Ys{g) are bounded with 
respect to the weak metric (3) since the functions Aj., Bk, k G Z^ U {0}, 
are bounded on T^, and y{- , ■ ) is bounded on [t,T] x T^. Moreover a{g) is 
C°°-smooth and Yg is at least C^-smooth on G°. D 

One can also consider (11) as an SDE with values in the Hilbert space 

i/"(T2,M2)_ 

Theorem 3. There exists a unique strong solution Z^ to the iir"(T^,R^)- 
valued SDE (11) on [t,T], with the initial condition Z^^ = e where e is the 
identity of G^ . This solution coincides with the solution to SDE (11) on G" 
orG"". 

Proof. By Theorem 1, SDE (11) on G^ has a unique strong solution Z*'*^ on 
[t,T]. Let us prove that the solution Z*'^ to (11) solves this SDE considered 
as an SDE in if"(T^,M^). Consider the identical imbedding ly : G" — >• 
if"(T^, M^), g h^ g. Applying Ito's formula to ly, and taking into account that 
AkigWig) = ^Ak^°g = ^kig) and that Ak{g)Ak{g)tvig) = Ak{g)Ak{g) = 0, 
we obtain that the solution Z*''^ to (11) on G° solves the iJ°(T^,M^)-valued 
SDE (11). Note that by the uniqueness theorem for SDEs in Hilbert spaces, 
SDE (11) can have only one solution in L2(T^, M^). This proves the uniqueness 
of its solution in //"(T^R^) as well. Thus the solutions to (11) on G", G^, 
and in //"(T^M^) coincide. D 

Let us find the representations of SDE (11) in normal coordinates on G" 
and Gy . First, we prove the following lemma. 



11 



Lemma 9. The following equality holds: 

/s rs 

i.e. instead of the ltd stochastic integral in (11) we can write the Stratonovich 
stochastic integral J^ a{Z^'^) o dWr- 

Proof. We have: 

a{Z'/) o dWr = a{Z'/)dWr + Y^ dAk{Z','')df3i{r) + dBk{Z'/)df3i{r). 

fcgZ+U{0},|fcK7V 

Hence, we have to prove that rfAfc(Z*'")d/3^(r) = and dBk{Zl'^)d(3^{r) = 0. 
For simphcity of notation we use the notation A^ for both of the vector fields 
Ak and Bk and the notation A^ for Ak and Bk, k G Zj' U {0}. Also, we use 
the notation f3^{s) for the Brownian motions {f3^{s), /3^(s)};.gg+Li|Q| i^i^j^- We 
obtain: 

d{A, o zf ) = J2 Mzln {A. o z'/) o d/3,(s) + y/'^(A, o z'/)dt. 

7 

This implies 

d{A, o Z'/) ■ df3, = A,{Z'/) {A, o Zl'')ds = 

which holds by the identity [k, V) cos(A; ■ 6) = {k, V) sm{k ■ 9) = or by 
differentiating of constant vector fields. D 

Let Zl = {Z*'^'^, ^s'^'^}fcgz+u{o} t>e the vector of local coordinates of the 
solution Z*'^ to (11) on G", i.e. the vector of normal coordinates provided 
by the exponential map eicp : T^G^ — )■ G". Let Ue be the canonical chart of 
the map eicp. 

Theorem 4 (SDE (11) in local coordinates). Let 

r = inf{sG[t,T]:Zf ^f/J. (12) 

On the interval [t,T], SDE (11) has the following representation in local 
coordinates: 



Ztr = P yr^'^dr + 4e {Pt{s A r) - /3,^(t)), 

Z*xf = ^'"' Y^-^'^dr + 6,e (/3f (s A r) - /3f (t)). 
w/iere 5k = I if \k\ ^ iV, anc? 6k = if \k\ > N. 

12 



(13) 



Proof. Let g = {g^^,g''^}k£z+u{o} ^^ local coordinates in the neighborhood 
Ue provided by the map eicp. Let / G C°°(G"), and let / : T^G" — )■ M be such 
that f = f oexp. Since eicp is a C°°-map (see [10]), then / e C°°{Uo), where 
Uo = eip-if/e. Note that g^f{g) = Akig)f{g) and g^fig) = B,,{g)f{g). 
By Ito's formula, we obtain: 

f{zitr) - /(e) = KzZ) - m 

"SAt Qf psAt 



/sAt Of psAt a/' 

* E |4f(2;)n'"+/ 'E 4^(z,')n'"o<(r) 

fcez+u{o} fcez+u{o} 

fcGZ+U{0} fcGZju{0} 

= r^dr V (F/^'=^A,(Z*''^)/(Z:'^) + F/^'=^i?,(Z*'^)/(Z;'^)) 



* fcGZju{0} 



sAt 



e 5^ 5, (A,(Z*'^)/(Z*''^) o dP^ir) + B,{Zlnf{Zln ° ^/^f (r))- 



* fcGZ+U{0} 

Using representations (10) and (6) we obtain: 

sAt psAt 



f{Zltr)-f{e)= r \r{Z'/)f{Z'/)dr+ f \a{Zl'^)f{Z'/)odWr. 
Jt Jt 



H Jt 

This shows that the process 



exp{ Y. Zlt^A, + Zlt^B,] 



fcGZ+U{0} 

solves SDE (11) on the interval [t, r]. D 

Let 

7t _ r yt;kA r^t-kB yt;kA yt;kB ytfiA ^t;OB\ 

^S — i^S 5^S y^S 5^S 5^S 5^S I k&Z+ 

be the vector of local coordinates of the solution Z*'"^ to (11) on G"^, i.e. the 
vector of normal coordinates provided by the exponential map exp : TgG" — )■ 
G". Further let Ue be the canonical chart of the map exp. 



13 



Theorem 5. Let 

f = inf{se[t,T]:Zf ^f/e}. 

Then, a.s. f = r, where the stopping time r is defined by (12), and on [t,r], 
Zt-,kA ^ ^t-MA^ zf^ ^ ^t;feB^ kellU {0}, Zf^ = Zf^ = 0,ke Z+, a.s. 

Proof. Let us introduce additional local coordinates g''-^,g^^, k G Z^, and 
perforin the same computation as in the proof of Theorem 4. We have to 
take into account that Y^-^ = Y^^ = 0, /c G Z^, and that the components of 
the Brownian motion are non-zero only along divergence-free and constant 
vector fields. We obtain that the coordinate process Z* verifies SDEs (13) 
and the equations Z^^''^ = Z^^^^ = 0, keZJ. D 

4.2 The backward SDE and the solution of the FBSDEs 

We have the following result: 

Theorem 6. Let Yg be the right-invariant vector field generated by the solu- 
tion y{s, ■) to the backward Navier-Stokes equations (9). Further let Zl'^ be 
the solution to SDE (11) on G^. Then there exists an e > such that the 
triple of stochastic processes 

Z*'^ r/'^ = tiZln. X'f = ea{ZY)%{Z'f) 
solves FBSDEs (5) on the interval [t,T]. 
Remark 3. The expression a{Zl'^)Ys{Z]:'^) means the following: 



a(zf )y;(zf ) = Y. MzTmzT) ®4 + B,{zlnuzln ® e 

k€Z+U{0},\k\fiN 



where Ys{-) is regarded as a function G" — ?■ iJ"(T^,M^), and Ak{g)Ys{g) 
means differentiation of Yg : G" — )■ iJ"(T^,M^) along the vector field A^ at 
the point g G G" • Let 7^ be the geodesic in G" such that 70 = e and Jq = Ak. 
We obtain: 

A,ig)ti9m = -^til^og)ie)\^^, = R^-^yis,^Sk=0 

= R,VA,y{s,e) = VAM9){e). (i4) 
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Thus, 

XT = e Y. \^A,y{s. ■)®et + Vs,y{s. ■ ) ® ef ] o Zf , (15) 

A:eZ+U{0},|fe|s£Af 

and the stochastic integral in (7) can be represented as 

i-T 

X'/dWr 

= e E / ^A,y(r, ■ ) o Z'/d(5i{r) + I VsMr, ■ ) ° Z'/d(5i{r). 

fcGZ+U{0},|fc|^7V 

In particular, if A^ = 0, 

X^dWr = e (^j^ J-2/(r, ■ ) o Zf d/3o^(r) + ^ J-y(r, ■ ) o Zf d/3o^(r) 

A result similar to Lemma 10 below was obtained in [6]. 

Lemma 10 (The Laplacian of a right-invariant vector field). Let V be the 
right-invariant vector field on G"' generated by an H°'~^^ -vector field V on 

T^. Further let e > be such that % ( 1 + | Sfcez+ \k\^N Ifcp^ ) ~ ^- Then for 

all g G G^, 

Y E {^A,VA,+VB,VB,)V{g) = uAVog. (16) 

A;eZ^U{0}, 

Here a is an integer which is not necessary equal to a. 

Proof. By the right-invariance of the vector fields Va^. V^^V^ and Vb^. V^^,!^ 
(Lemma 7), it suffices to show (16) for g = e. We observe that 

(fc, V) cos(A; ■ 6) = {k, V) sin(A; ■ ^) = 0. 

Then, for A; e Z+, ^ G T^, 

^A,\^A,viem = —L_cosik.e)Ck,v)[cosik.e)Ck,v)vie)] 



lU2a+2 
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^ cos{k-e)\k,vyv{e). 



Similarly, VB,VB,t>(e)(^) = ^^^sm{k ■ eY{k,VfV{e). Hence, for each 

k e Z+, 

(Va,Va, + VB.VijJ V(e)(^) = ^^(fc,V)V(^). (17) 

Note that for each k G Z^, either A; or —A; is in Z^, and 

(^,V)2 + (A;,V)2 = |fc|2A. 

Summation over /c G Z^, |/c| ^ A^, in (17), and coupling the terms numbered 
by k and k (or —k) gives: 

^ (VA,VA, + VB,VBjV^(e)(^) = i Y. ^^^(^)- 

Note that (VaoV^o + V 3,^ Bo)V {e){d) = AV(9). Finally, we obtain: 

E (VA,VA,+VB,VBjV^(e)(^)=(l + i 5^ ^)^^(^)- 

fcGZ^UJO}, fcGZ^,|fc|^Ar 

\k\siN 

The lemma is proved. D 

Corollary 1. Let the function </? : T^ -)■ M^ fee C'^ -smooth. Further let 
Ak{g)[ipog] and Bk{g)[(pog], fc G Z^, mean the differentiation of the function 
G" — )• i^2(Tr^, I^^)) g ^-^ '^ ° g along A^ and resp. B^. Then for all g G G", 



e^ 



E {Ak{g)Ak{g) + Bk{g)Bk{g))[ipog] = v/\^og. (18) 



2+U{0}, 



Proof. The computation that we made in (14) but applied to if o g implies 
that 



Ak{g)[y^og]= \^-L_cos{k-9)Ck,VM9) 



'9- 



Similarly, we compute Bk{g)[ip o g]. Now we just have to repeat the proof of 
Lemma 10 to come to (18). D 
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Lemma 11. Let ^r, r e[t,T],t e [0,T), be an H°'{T'^,R'^) -valued stochastic 
process whose trajectories are integrable, and let (pT be an H°'{T'^,M.'^) -valued 
random element so that both $r and(f)T possess finite expectations. Then there 
exists an J^g-adapted if"(T^, M^) xC[E, iJ°(T^, M^)) -valued pair of stochastic 
processes {Ys,Xs) solving the BSDE 

Ys = (t)T+ I $r dr- j Xr dWr (19) 

J s J s 

on [t,T]. The Yg-part of the solution has the representation 

T 



r, =E[0T+ / $rc/r|j-,], (20) 

J s 

and therefore is unique. The Xg-part of the solution is unique with respect to 
the norm ||^s|| = J^ ll^sll£(i<;,_H'™(T2,i82)') ds. 

The proof of the lemma uses some ideas from [20]. 

Proof. Representation (20) follows from (19). Let us extend the process Yg 
to the entire interval [0, T] by setting Ys = Yt for s G [0, t], and note that the 
extended process Yg is a solution of the SDE 

Ys = (t)T+ I \t,T] $r dr- I XrdWr 
J s J s 

on [0,T]. Let X, G C{E,H'^{T'^,M.^)), s G [0,T], be such that 

^[<Pt+ [ h,T] $r dr-Yo\ J-,] = f Xr dWr. (21) 

Jo Jo 

The process Xs exists by the martingale representation theorem. Indeed, on 
the right-hand side of (21) we have a Hilbert space valued martingale. 

By Theorem 6.6 of [12], each component of the i7"(T^,M^)-valued 
martingale on the right-hand side of (21) can be represented as a sum 
of real-valued stochastic integrals with respect to the Brownian motions 
{Pki^)^ Pk{^)}k&'L+(j{o\\kKN- Hence, there exist J^^-adapted stochastic pro- 



cesses {X,^^,Xi'^}fcez+u{o},|feK7v such that 

E[0T+/" \t,T]^rdr-Y^\Fs] = V r X^''d(3^ir)+ fx,^^ rf/3f (r). 

JO , ^1 r„i JO JO 



fcGZ+U{0}, 

\k\<:N 
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Let the process Xg be defined by (8) via the processes X^^ and X^^, k G 
2,2 U {0}, 1^1 ^ N. Ito's isometry shows that E J^ Il-^r||£(£;/^c,(']r2 r2)) < oo. 

Note that for all s G [0,t], ^ X^dWr = f^X^dWr- This shows that X, = 
for almost all w G ^2 and almost all s G [0,t], and therefore can be chosen 
equal to zero on [0,t]. Thus, (21) takes the form: 

E [0T + / ^rdr-Yt\J^s]= f Xr dWr. (22) 

It is easy to verify that the pair {Ys,Xs) defined by (20) and (22) solves 
BSDE (19). To prove the uniqueness, note that any J^,-adapted solution to 
(19) takes the form (20), (22). Moreover, if the processes Xg and X'^ satisfy 
(22), then 



y^a ^o r/c uarir2 m2\\ UT 



C{E,H°'(T^ 



T 



{Xs - X',) dWr 



2 

= 0. 



n 



Proof of Theorem 6. Let us consider BSDE (7) as an L2(T^,M^)-valued 
SDK, and Yg as a function G" — )■ L2(T^,M^). Since for each s G [t,T], 
vis, ■ ) G H'^+^T^, M2) and a > 2 by assumption, then y, : G ° -> ^2(1^, R^) 
is at least C^-smooth. Equations (2) show that the function dsy{ ■ , ■ ) : 
[t,T] — )■ L2(T^,M^) is continuous since Vp, Ay, and {y,'Vy) are continu- 
ous functions [t,T] — )■ L2(T^,M^) by Assumption 1. Taking into account 
that the diffeomorphisms of Gy are volume-preserving, we conclude that for 
each fixed g G G", dsYs{g) : [t,T] — )■ L2(T^,M^) is a continuous function. 
Hence, Y, : [t,T] x G^ ^ L2{T\R'^) is C^-smooth in s G jt,T] and C^- 
smooth in (yf G G". Ito's formula is therefore applicable to Ys{Zl'^). Below 
we use the fact that Z*'*^ is a solution to forward SDE (11) and the identity 
^(Z*''^) = q'^ °^l'^- For the latter derivative we substitute the right-hand 
side of the first equation of (2). The notation X {g)\Ys{g)] (sometimes without 
square brackets) means differentiation of the function Yg : G° — > L2(T^,]R^) 
along the right-invariant vector field X on G" at the point g G G" . The 
same argument as in Remark 3 implies that X{g)[Ys{g)] = VxYs{g). Taking 



into account this argument, we obtain: 
t{Z'f)-h{Ztf) = - I dX{Z'/)dr- I dr%{Z'/)\%{Zln] 

J s J s 

- I ea{Z'/)Yr{Z'/)dWr. (23) 

J s 

Note that 

UzlnitiZtn] = [(2/(r, ■),V)y(r, .)]oZt'^. 
Also, let us observe that 

^ E [MKnAkiztntiztn + B^izinB^iztntiztn] 

feeZju{0},|fe|^Ar 

= u[Ay{s, ■ )] o Z*''^ 

where the latter equality holds by Lemma 10, and e > is chosen so that 

t(i + lEfcGZ+,|fcKiv^) = '^^ Note that the terms Va,Va,V;.(^*'^) and 

V B^^ Bk'^riZl''^) are elements of TG""^, and therefore are well defined in 
L2(T2,M2J_ Continuing (23), we obtain: 

YMn - Hz'tI 

rp 

= j dr [V^(r, Z'/) + [(y(r, ■ ), V)y{r, ■ )] o Z*'^ + u[Ay{r, ■ )] o Z*'^" 
- / [(y(r, ■),V)y(r, OloZ^'^rfr- f u[Ay{r, ■)]o Z'/ dr 

J s J s 

- f ea{Z'/)YriZ'/)dWr= f V{r,Z'/)dr- I ea{Z'/)%{Z'/) dWr- 

J s J s J s 

(24) 
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Thus the pair of stochastic processes {Ys{Zl''^),ea{Zl''^)Ys{Zl''^)) is a solu- 
tion to BSDE (7) in L2(T^,M^). It is J-'^-adapted since Z*'^ is J-'^-adapted. 
By Lemma 11, we know that there exists a unique J-'^-adapted solution 
(F/'^X*•^) to (7) in if"(T2,M2). Clearly, (y/•^X*'^) is also a unique J",- 
adapted solution to (7) in L2{T^,M?). Hence, Y^^'^ = Ys{Zl'^) and /f ||X*'^ - 
ea{Zl''^)Ys{Zl''^))\\'^,^ jja,j2^2\)ds = 0, and therefore the pair of stochas- 
tic processes (1^(Z*'^), e(j(Z*''^)l^(Z*''^)) is a unique J-'^-adapted solution to 
BSDE (7) in //"(T^, R^). The theorem is proved. D 

5. Some identities involving the Navier— Stokes solution 

The backward SDE allows us to obtain the representation below for the 
Navier-Stokes solution. Also, it easily implies the well-known energy identity 
for the Navier-Stokes equations. 

5.1 Representation of the Navier-Stokes solution 

Theorem 7. Lett e [0,T], and let Z^ he the solution to SDE (11) on [t,T] 
with the initial condition Z^'^ = e. Then the following representation holds 
for the solution y(t, ■ ) to (9). 

y{t, ■)=E [h{Z'/) + J Vp{s, ■ ) o Z'/ds 

Proof Note that Yt{Zt'^) = vit, ■), and E[J^ X^^'^dWr] = 0. Taking the ex- 
pectation from the both parts of (7) at time s = t we obtain the above 
representation. D 

5.2 A simple derivation of the energy identity 

Ito's formula applied to the squared L2(T^,]R^)-norm of Y^''^ gives: 

WY^m = \\h{Z'/)\\l 
+ 2 f {Y:'\V{Z'/))L,dr-2 f {Y:'\X'/dWr)L,- f \\X'/\\ldr. (25) 

Js J s J s 
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Using representation (15) for the process X*'*^ we obtain: 

fc6Z+U{0},|fc|<Af 



Ilz 



'' l\\L2 



e 



>-)llL 



^'(l + ^ E ^)llVy(^,-)llL = 2H|Vy(s,.)ll^ 



2 f- |A;|^ 

A:GZ^,|fc|^Ar 

Taking the expectation in (25) and using the volume-preserving property 
of Z*'^, we obtain: 

lb(^, ■)llL + 2^rilv^(^'-)lliA = ||/^|ii^. 



6. Constructing the solution to the Navier— Stokes equa- 
tions from a solution to the FBSDEs 

Let us prove now a result which is, in some sense, a converse of Theorem 6. 
In this section we consider (5) as a system of forward and backward SDEs in 
the Hilbert space ilf"(T^,R^), where a ^ 3. As before, let V^(s,Z*'^) denote 
Vp(s, ■ ) o Z*'*^, and let J^s denote the filtration a{Wr,r G [0, s]}. 

Theorem 8. Assume, for an H^'^^-smooth function p{s, ■ ), s G [0,T], and 
for any t G (0, T), the existence of an J^g-adapted solution [Z]:^, y^*'*^, X*'^) to 
(5) on [t, T] such that the processes Z^ andY^''^ have a.s. continuous trajec- 
tories and such that Z^'^ take values in G^. Then there exists Tq > such 
that for all T < Tq there exists a deterministic function y{s, ■) G T^G^ on 
[0,T], such that a.s. on [t,T] the relation Y^''^ = y{s, ■ )oZl''^ holds. Moreover, 
the pair of functions {y,p) solves the backward Navier-Stokes equations (9) 
on [0,T]. 

Lemmas 12-18 below are the steps in the proof of Theorem 8. 
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Lemma 12. For all t G [0,T) and for any J^t-f^Gosurahle G^-valued random 
variable C,, the triple of stochastic processes 

(Z'/, Y'/, X'/) = {ZY o e, n*''^ o e, XY o (26) 

is J^s-O'dapted and solves the FBSDEs 



Z'J^ = e + /; y/ dr + /; a{Z'/) dW, 



(27) 



on the interval [t,T] in the space H°'{T^ 



Proof. Let us apply the operator R^ of the right translation to the both sides 
of FBSDEs (5). We only have to prove that we are allowed to write R^ under 
the signs of both stochastic integrals in (5). Let us prove that it is true for 
an J^i-measurable stepwise function ^ = ^'^iQi'^A,, where gi G G" and the 
sets Ai are J^r measurable. Indeed, let s and S be such that t ^ s < S ^ T, 
and let $r be an J>-adapted stochastically integrable process. We obtain: 

pS oo oo „s °° pS 

/ <^r dWr °^9dA = Yl ^^» / '^r°gidWr = ^ Ia, $r O Qi dWr 

i-S °° 

= / <^roJ29dA, dWr. 
Js ■_^ 



i=l 



Next, we find a sequence of J-i-measurable stepwise functions converging to 
^ in the space of continuous functions C(T^,]R^). This is possible due to 
the separability of C(T^,M^). Indeed, let us consider a countable number of 
disjoint Borel sets O^ covering C(T^, M^), and such that their diameter in the 
norm of C(T2,M2) is smaller than i. Let A^ = ^"^0,") and ^f G O"} n G^. 
Define ^„ = YllLi 9i ^A"- Then it holds that for all u & Q, ||^— ^n||c(T2,M2) < -. 
Let /($) and /($ o ^) denote /^ $r dWr and resp. /^ ^r°^ dWr- We have to 
prove that a.s. /($) ^ S, = /($ o ^). For this it suffices to prove that 

hm E||/($) o e„ - /($) o e||i^(Tr.,j,.) = 0, (28) 

hm E||/(<|. oe„) - mo^)\\l , , = 0. (29) 

Due to the volume-preserving property of ^ and ^n, \\H^) ° C«IIl2(t2r2) ~ 
ll^(*) °'^IIL(T2,R2) = II^(*)IIL(T2,R2)- Hence, by Lebesgue's theorem, in (28) 
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we can pass to the limit under the expectation sign. Relation (28) holds then 
by the continuity of /($) in 6* G T^. To prove (29) we observe that by Ito's 
isometry, the limit in (29) equals to lim„_^ooIE J^ ||$j. o^„ — $^0^11^ ^ ^2\dr. 
The same argument that we used to prove (28) implies that we can pass to 
the limit under the expectation and the integral signs. Relation (29) follows 
from the continuity of $r in 6' G T^. 

Hence, (Z*'"^ o ^,1^^*'"^ o ^jX*'"^ o ^) is a solution to (27). This solution is 
clearly J^^-adapted. D 

Lemmas 13-17 below use some ideas and constructions from [9]. 

Lemma 13. The map [0,T] x T^ ^ R'^, {t,e) h-> F/'"(^) is determmistic. 

Proof. Let us extend the solution (Z*'^, Y'^*''^,X*'^) to the interval [0,t] by 
setting Z*'*^ = e, F/'"^ = Y^ '^ , Xl'^ = for all s G [0,t]. The extended process 
solves the problem: 



'l'^ = e + /; I[t,T] (r)F,*'^ dr + /^ I[,^] (r)a(Z*.^) dWr 
n*'^ = hiZ'/) + /Jl[,,r](r)V^(r, Z^'-) dr - jJx'^'-dWr. 



The random vector Y^ '^ is J-Q-measurable, and hence is deterministic by Blu- 
menthal's zero-one law. Since Yf- '^ = Yq''^, the result follows. D 

Lemma 14. There exists a constant Tq > such that for T < Tq the function 
[0,T] -> if2(T2,R2), t ^ yI'^ is continuous. 

Proof Let (Z*'^ F,*'^ X*'^) and (Z*'^ y/'^ X*'"^) be solutions to (27) which 
start at the identity e at times t and resp. t', and let t < t' . These solutions 
can be regarded as solutions of (30) if we extend them to the entire interval 
[0,T] as it was described in Lemma 13. The application of Ito's formula to 
||y^*'^||^ .,jj,2 JJ2-) and the backward SDE of (27) imply that the expectation 
^II^s^'^^IIl (t2 MS) i^ bounded. The forward SDE of (30), Gronwall's lemma, 
and usual stochastic integral estimates imply that there exists a constant 
Ki > such that 



E||Z*'^ - Zl'''\\l^(j2^^2) < Ki 



'-Jo 



Let us apply Ito's formula to 11^/'*^ — ^s*''^|li2(Tr2 k2) when using the backward 
SDE of (30). Again, Gronwall's lemma, usual stochastic integral estimates 
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and the above estimate for E||Z*'^ — ^^'^||^ (^2 ir2) imply that there exists a 



constant K2 > such that 

i-T 



E||F^^-K*'1lL.T2,.2^<ir, 






We take Tq smaller than -^. Then there exists a constant i^ > such that 

sup E||n*'^ - r/'^||i,(Tr2,K^) < Kit' - t). (31) 

se[o,T] 

Evaluating the right-hand side at the point s = t, and taking into account 
that Y^ '^ = Y^, ''^ we obtain that 

ll^t' ~^t'' IIl2(t2,ir2) < -^"^(^ ~^)- (32) 

Differentiating (30) with respect to 9 we obtain the following system of for- 
ward and backward SDEs: 

'VZl'^ = / + /; I[t,T](r)Vy,*'^ dr + /; I[,,T](r)Va(Z*'^) VZ*'^ dWr 
Vy,*'5 = V/i(Z^")VZ^" + /Jl[i,T](r)Vy(r, Z*'5)VZ*'^rfr 

Again, standard estimates imply the boundedness of E||VZ*''^||^_^™2 ]g2) and 
E||Vy/''^||^^/.]j.2 J52)- The same argument that we used to obtain (32) as well 
as the estimate for the sup^g^o.r] JE||2'*''^ — Zl'^\W,j2 ^2), which easily follows 
from (31), and the forward SDK imply that there exists a constant L > 
such that for all t and t' from the interval [0, T], 

\\VY^' — VF^,' ||l2(t2,r2) < I/|t —t\. (33) 

Differentiating (30) the second time and using the same argument once again 
we obtain that there exist a constant M > such that for all t and t' 
belonging to [0,T], 

II VVF/'^ - VVF;'^||i^(Tr2,i,2) <M\t'-t\. (34) 

Now (32), (33), and (34) imply the continuity of the map t i— )■ Y^ '^ with 
respect to the if ^^T^, M^). topology. D 
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Everywhere below we assume that T < Tq where Tq is the constant defined 
in Lemma 14. 

Lemma 15. For every t G [0,T) and for every Tt-measurable random variable 
^, the solution (Z*'^, F^*'^,X*'^) to (27) is unique on [t,T]. 

Proof. Let us assume that there exists another solution {Zl'^,Y^'^,Xl'^) to 
(27) on [t,T]. The same argument as in the proof of Lemma 14 implies the 
uniqueness of solution to (27). Specifically, the argument that we applied to 
the pair of solutions (Z*'% y/'% X*'^) and (Z*''^, F/''^,X*'''^) has to be applied 
to (Z*'^,F/'^,X*'^) and (Z*''^, Y'^*'^,X*'^), and it has to be taken into account 
that t = t'. D 

Lemma 16. Let the function y : [0,T] xT"^ ^M."^ be defined by the formula: 

y{t,9) = Yl'%9). (35) 

Then, for every t G [0,T], y(t, ■ ) is H"'-smooth, and a.s. 

Y^'^ = yiu, ■ ) o Zf . (36) 

Proof. Note that (26) implies that if ^ is J-'t-measurable then 

Y!^^ = y{t,-)o^. (37) 

Further, for each fixed u G [t, T], (Z*'*^, Y^''^, X*'^) is a solution of the following 
problem on [m, T]: 

'Zl'^ = Z'^ + £ Yl^'dr + /J a{Zl^^)dW, 
Yl^' = h{Z'/) + J^V{r, Zl:')dr - /Jx*'^diy,. 

By uniqueness of solution, it holds that Y^'*^ = Y^'^^ a.s. on [m, T\. Next, by 
(37), we obtain that Y^'^^ = y{u, ■ )oZlf. This implies that there exists a set 
Vtu (which depends on u) of full P-measure such that (36) holds everywhere 
on VLu- Clearly, one can find a set ^q, P(f2(Q) = 1, such that (36) holds on 
VLq for all rational u G [t,T]. But the trajectories of Z*'^ and Y^''^ are a.s. 
continuous. Furthermore, Lemma 14 implies the continuity of y{t, ■ ) in t with 
respect to (at least) the L2(T^,M^)-topology. Therefore, (36) holds a.s. with 
respect to the L2(T^,M^)-topology. Since both sides of (36) are continuous 
in ^ G T^ it also holds a.s. for all ^ G T^. D 
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Lemma 17. The function y defined by formula (35) is C^ -smooth in t G [0,T]. 
Proof. Let 6 > 0. We obtain: 

yit + 6,-)-yit,-) = Yl^, 






■t+5 



t+S'-^ t+5 ^ t 



Let Yg be the right-invariant vector field on G"' generated by y{s, ■ ). Lemma 
16 imphes that a.s. 

Thus we obtain that a.s. 

y{t + 6,-)- y{t, ■ ) = (l>+,(e) - Y.^siZl':^)) + (Y^^s " yI'I- 

We use the backward SDE for the second difference and apply Ito's formula 
to the first difference when considering Yf^s as a C^-smooth function G" — )■ 

L2(T2,M2)_ We obtain: 



Yt+siZlts) - Yt+sie) 
-t+s 

drt'^^iZ'/)[Y,+siZ'/)] 

ft+S 



t+5 



ea{Z'/)Yt+s{Z'/)dW, 



dr J2 [MztnMztn+MztnBk{zin]yt+s{zin- 



fcez+u{o} 
The same argument as in Theorem 6 implies: 



Yf.x(Z,\.) - YtA-sie 



■t+s\^t+5 



t+5\ 



t+5 



drVyfr,-)yit + 5, ■)o Z. 



t,e 
r 



/t+5 j-t+5 

druAy{t + S,-)o ZY + / e o{ZY) Yt+siZ'/) dW,. 



'•t+5 



>'t+5 



Further we have: 
Finally we obtain that 



dWr. 



l{y{t + 5,.)-y{t,.))=--^¥. 



t+5 



dr[{y{r, ■),V)y{t + 6, ■) 
+ ijAy{t^6,-)^Vp{r,-)]oZ'/. (38) 
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Note that Z*'"^, Vp(r, ■ ), and {y{r, ■),'V) y{t + 5, ■ ) o Zl'^ are continuous in 
r a.s. with respect to the L2(T^,M^)-topology. By Lemma 14, V|/(t, ■) and 
A i/(t, ■ ) are continuous in t with respect to the L2(T^, ]R^)-topology. Formula 
(38) and the fact that Z^'^ = e imply that in the L2(T^,M^)-topology 

%(t, ■) = -[V,(t,.)2/(t, ■) + ^Ay(t, O + V^t, ■)]• (39) 

Since the right-hand side of (39) is an iJ°~^-map, so is the left-hand side. 
This implies that dty{t, ■ ) is continuous in 6* G T^. Relation (39) is obtained 
so far for the right derivative of y{t, 9) with respect to t. Note that the right- 
hand side of (39) is continuous in t which implies that the right derivative 
dty{t, 9) is continuous in t on [0, T). Hence, it is uniformly continuous on every 
compact subinterval of [0, T). This implies the existence of the left derivative 
oiy(t, 9) in t, and therefore, the existence of the continuous derivative dtyit, 9) 
everywhere on [0,T]. D 

Lemma 18. For every t G [0, T], the function y{t, ■ ) : T^ — t- R^ is divergence- 
free. Moreover, the pair {y,p) verifies the backward Navier-Stokes equations. 

Proof. Fix a t > 0, and consider the TgC"- valued curve 7^ = E [exp~^Z/^], 
C ^ t, in a neigborhood of the origin of T^G^. The forward SDK of (27) can 
be represented as an SDE on G": 

'dZl'' = exp{f;(Z*'^) ds + aiZ'f) dWs}, 
Zt = e, 

where 1^ is the right-invariant vector field on G" generated by y(s, ■ ). This 
implies that 

d 



dC^^ 






and therefore y{t, ■ ) G TgG" . Next, the backward SDE of (27) implies that 
Yrp'"^ = h{Zr^^). This and relation (36) imply that y(T, ■) = h. Since we 
already obtained (39) in Lemma 17 the proof of the lemma is now complete. 

D 

7. The backward SDE as an SDE on a tangent bundle 

Let (Z*•^y,*•^X*''^) be a solution to FBSDEs (5). We will show that the 
backward SDE can be represented as an SDE on the tangent bundle TG" as 
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well as an SDE on TG". We will construct a backward SDE in the Dalecky- 
Belopolskaya form (see [4]) and show that the process F/'^ is its unique 
solution. 

7.1 The representation of the backward SDE on TG" 

Let y{s, ■ ), s G [t,T], be the solution to the backward Navier-Stokes equa- 
tions (9). Let Yg be the right-invariant vector field on Gy generated by y{s, ■ ). 
The connection map on the manifold Gy generates the connection map on 
the manifold TGy as it was shown in [4], p. 58 (see also [11]). As before, 
we consider the Levi-Civita connection of the weak Riemannian metric (3) 
on Gy. Let exp denote the exponential map of the generated connection on 
TGy. More precisely, exp is given as follows: 

"^ _ /^7a(l) 



where ( ,,! ] is the geodesic curve on TG" with the initial data 7' (0) = a, 

?7^(0) = /3, 7a (0) = X, 77/3(0) = a. Let the vector fields A^ and Bj^ be the 
horizontal lifts of A^ and B^ onto TTG^. Further let (9,y/ be the vertical 



V ' 



lift of dgYs onto TTG". Let us consider the backward SDE on TGy-. 

dY'^^' = expyt.e^dsY^\Y^'^')ds + S{Y^'^'')ds 

+ e Yl [A^{Y:n®et + B^{Y:n^et]dW,], (40) 

fc6Z+U{0},|fc|^Ar 

where S is the geodesic spray of the Levi-Civita connection of the weak 
Riemannian metric (3) on Gy (see [15] or [16]), and Zl'^, s G [t, T], is the 
solution to (11) on Gy with the initial condition Z^''^ = e. 

Theorem 9. There exists a solution to (40) on [t,T]. Moreover, if 
dsy{s, ■) G if°(T^,]R^), then this solution is unique and coincides with the 
Y^'^-part of the unique J^s-O'dapted solution (y/'^,X*''^) to (7). 

Proof. From the proof of Theorem 6 we know that the pair of stochastic pro- 
cesses {Ys{.Zl'^) , e a{Zl'^)Ys{Zl'^)) is the unique J%j-adapted solution to (7) in 

28 



if"(T^, M^). Let us prove that Ys{Zl'^) is a strong solution to (40). First we de- 
scribe a system of local coordinates {g^^ , X^^ , g^^ , X^^) ^^■^+yJsn^ in a neigh- 
borhood UegxTeG^ of the point X{g) G TG^ where f/g C G" is the canonical 
chart. The vector g = {g^^, g^^)kez+uio} ^^ ^^^ vector of normal coordinates 
in the neighborhood Ueg, g G G" . The vector X = {X^^, -^^^)fcez+u{o| repre- 
sents the coordinates of the decomposition of the vector X{g) G TG" in the 
basis {A,,5fc},,^+^{oj: X{g) = j:,^^^^^,^{X'^A,{g) + X''^B,{g)). Let / be 

a smooth function on TG", and let f{X,g) = f{X{g)), where X{g) G TG". 
Let r be the exit time of the process Z*'*^ from the neighborhood UgZl''^. 
We will compute the difference f(Y^''^) — f(Y^'^) using Ito's formula. Let 
{Zr, Yr) = {Z^^, Z^^, Y^^, ^r^)k&+vj{o} ^^ ^^^ vector of local coordinates of 
the process Yr{Z^''^) on [s,t]. Using SDE (40), we obtain: 

^,df{Yr,Zr) , ,,,kBydm,Zr) 



m'n-fiY:n = - E / i^'"')' ^.^iT + (^^ 



fcez^u{o} 



dYj'^ ■ V r ; QYJ 



kB 



' r 



QZkA ^r g^kB ^ 2''\d{Z^^y d{Z^^)V ^ ^ ""' 



Zjr- 



dr 



d f{Yr,Zr) , dfiY^Zr) ^ ^^ 



^ E / ^^^^9^^e^+-A^^^e^.iy. (41) 

fcgZ+U{0},|fc|s£Af * 



where (5^ = 1 if |/i;| ^ N , and 5fc = otherwise. Since / is a smooth function 
on TG^, all its restrictions to the tangent spaces of G^ are smooth. Hence, 
one can talk about derivatives of / restricted to a tangent space along the 
vectors of this tangent space. Namely, the following relation holds: 

^£^ = fitizt-))Mzn 

Note that the differentiation of / with respect to Z^"^ and Z!^^ can be re- 
garded as the differentiation of the composite function foYr along the vectors 
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Aj, and Bu. Namely, ^^^ = A^(Z*'^)[(/ o %){Zl'% This implies: 
f{Yln - mz'/)) = - r dr [drif o Yr)iZtn + tiztnif o %){Zln 

J s 

fcgZ^U{0},|fc|^7V 

fcgz+u{o},|fc|s;Af * 

(42) 

We extended the integration to the entire interval [s,T] since the local coor- 
dinates no longer appear under the integral signs. This is also possible since 
(41) holds also with respect to the local coordinates in the neighborhood 
UeZl''^ and a new exit time ri. The same argument can be repeated with 
respect to the local coordinates in the neighborhood UeZ!p^, etc. Let us con- 
sider now f o Yg as a time- dependent function of (7 G G". Applying Ito's 
formula to {f oYs){Zl'^) on the interval [s,T] and using SDE (11) on G", we 
obtain exactly the above identity. This proves that Y"/'*^ = Ys{Zl'^) is a strong 
solution to (40) on TG^. By results of [15], dsY^ is C^-smooth. Moreover S, 
A^ and Bj} , k G Z^, are C°°-smooth. Again, by results of [15], the solution 
of BSDE (40) on TG^ is unique. D 

7.2 The representation of the backward SDE on TG° 

Applying Proposition 1.3 (p. 146) of [4] (see also [16], p. 64) to the manifolds 
TGy and TG" and the identical imbedding ty : TG^ — )■ TG", we obtain 
that the process Zv(Kj(Z*''^)) = Ys{ZY) solves the following backward SDE 
on TG": 

dYl'^ = expy..{dsY;{Y:nds + S(n*'^)rf. 

+ e E K{Y:n®et + B^{Y,'n®ei]dW,], (43) 

Y^'' = h{Z'/) 

where S is the geodesic spray of the Levi-Civita connection of the weak 
Riemannian metric on G", dsY^ denotes the vertical lift of dsYg onto TTG°', 
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Aj^ and B^ denote the horizontal hfts of A^ and Bk onto TTG°', the process 
Zf^, s e [t, T], is the solution to (11) on G" with the initial condition Zj'^ = e. 
The exponential map exp on TTG" is defined similarly to the map exp on 
TTG^. Namely, the Levi-Civita connection of the weak Riemannian metric 
on G" generates a connection on TG". The latter gives rise to the exponential 
map exp on TTG" as it was described in Paragraph 7.1. We actually have 
obtained the following theorem. 

Theorem 10. Backward SDE (43) has a unique strong solution. Moreover, 
this solution coincides with the unique strong solution to BSDE (40) on TG", 
and with the Y^'^^-part of the unique J^g-o-dapted solution [Y^''^ , Xl''^) to (7). 

Proof. We have already shown that the process Ys{Zl'^) solves BSDE (43). 
The uniqueness of solution can be proved in exactly the same way as the 
uniqueness of solution to (40) on TG^ (see the proof of Theorem 9). D 

8. Appendix 

8.1 Geometry of the group of volume-preserving diffeomorphisms 
of the n-dimensional torus 

Let T" = 5^ X ... X S^ denote the n-dimensional torus. Let us describe 

^ v ' 

n 

a basis of the tangent space T^G^ of the group G^ of volume-preserving 
diffeomorphisms of T". We introduce the following notation: 

Z+ = {(A;i, A;2, . . . , A;„) e Z" : A;i > or A;i = ■ • • = ki_i = 0, ki> 0, 

i = 2,...,n}; 



k — (/ci, . . . , kn) G /i„ , \k\ 



.Y.^1 k-e = J2ki0i, 



i=l 



^ = (^i,...,^„)GT", V=(,|-,^,..., ^ 



dOi '802 dOn 

For every k G Z^, (A;^, . . . , k^"^) denotes an orthogonal system of vectors of 
length \k\ which is also orthogonal to k. Introduce the vector fields on T": 

^l = 7^cos(A;-e)F, Bl = -^^sin{k-e)k\ z = l,...,n-l, A; G Z+, 
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and the constant vector fields Aq, i = 1, . . . , ra, whose ^th coordinate is 1 and 
the other coordinates are 0. Let A\,B'l, i = 1, ... ,n — 1, k G Z+, denote 
the right-invariant vector fields on G" generated hj A]^, B^, i = 1, . . . ,n — 1, 
k G Z^, respectively, and let Aq = Aq, i = 1, . . . ,n, stand for constant vector 
fields on G^. The following lemma is an analog of Lemma 6. 

Lemma 19. The vectors Al{g), Bl{g), k G Z+, i = 1, . . . ,n-l, g e G", A^, 
i = 1, . . . ,n, form an orthogonal basis of the tangent space TgG^ with respect 
to both the weak and the strong inner products in TgG". In particular, the 
vectors A\, B\, k G Z+, i = 1, . . . ,n — 1, Aq, i = 1, . . . ,n, form an orthogonal 
basis of the tangent space Tf,G^. Moreover, the weak and the strong norms 
of the basis vectors are bounded by the same constant. 

The other lemmas of Section 2 hold in the n-dimensional case, with re- 
spect to the system Al, B\, k G Z+, i = 1, ... ,n — 1, Aq, i = 1, . . . ,n, 
without changes. The index a of the Sobolev space if ° is an integer bigger 
than f + 1. 

8.2 The Laplacian of a right-invariant vector field on G"(T") 

One of the most important steps in the proof of Theorems 6 and 8 is Lemma 
10, i.e. the computation of the Laplacian of a right-invariant vector field on 
G" with respect to the subsystem {^fc, -Bfc}fcgz+u{o},|fe|s£Af where N can be 
fixed arbitrary. Below we prove an n-dimensional analog of this lemma. 

Lemma 20. Let V be the right-invariant vector field on G"(T") generated by 
an H"''^'^ -vector field V on ¥"■. Further let e > be such that 

e^ /, n- 1 sr^ 1 \ 

1 + Yl Tr^)=^- 



2 V n ^^ l/tP" 



Then for all g G G°', 



n-l 



Proof. As it was mentioned in the proof of Lemma 7, it suffices to consider 
the case g = e. We observe that for alH = 1, . . . , n — 1, 

{k\ V) cos(A; -6) = - sm{k ■ e){k\ k) = 0. 
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Similarly, {k\ V) sin(A; ■ 9) = 0. Then, for keZ+,9 e T"", 

n—l ^ n—1 

5^ V^. V^. V(e)(^) = ^^^ 5^ cos(A: ■ e)(^\ V) [cos(A; ■ e)(^\ V)n^ 

i=l ' ' i=l 

1 n-1 _ 

= ^^^ ^^<k ■ ef ^{k\ V) V(^) = ^^^ cos(fc ■ en\k\'A - {k, vrwie) 

The latter equality holds by the identity Y^'^'^ik^V)^ + {k,Vf = \k\^A 
that follows, in turn, from the fact that the system {m' t|t}' "^ ~ 1, . . . , n — 1, 
forms an orthonormal basis of M". Similarly, 

Y, v^^s^{e)ie) = jj^, Mk ■ en\k\'A - {k, vfWie). 

Hence, for each A; G Z+, 

n— 1 ^ 

E(^^I^a; + ^Bi^BOVi^m = m^^d^l'^ - ik,Vf)V{e). (44) 
Further we have: 

fceZ+,|fc|s£Af keZ„,\k\iiN 



n 



= 5 E ij4tjE*.'9?+ E ^E **«<«. 

fceZ„,|fc|^Ar ' ' i=l keZn,\k\^N ' ' i^j 

where di = -^, and due to the factor | we perform the summation over all 
k G Z„. Clearly, the second sum is zero. To show this, we have to specify the 
way of summation. Let us collect in a group the terms kikjdidj attributed to 
those A; G Z,„ whose coordinates except the zth and the jth coincide, while 
the zth and the jth coordinates satisfy the following rules: they are obtained 
from ki and kj attributed to one of the vectors of the group by means of 
an arbitrary assignment of a sign. This operation specifies four vectors. The 
other four vectors are obtained from the first four vectors of the group by 
means of the permutation of the ith and the jth coordinates. In total, we 
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get eight vectors in the group. Clearly, the summands kikjdidj attributed to 
these vectors cancel each other. Let us compute the first sum. 



Z^ |^|2q+2 Z^^i^i Z^ [ Z^ \U2a+2^i 



df 



Note that 



E *H--= E *s4 E 1*1^- 



n 

kalj-n ,\k\=const A;SZn,|A;|=const fcgZ,i,|fe|=const 



This implies: 



V ^ V ^-2/52 = 1 V ^— A = - V — 

Z_^ |;^|2a+2 Z_^ J"^* „ Z^ |i,|2Q „ Z^ |^|: 

Together with (44) it gives: 






n-1 



E E(Va,v,.+v,. v,.)f (e)(^) = ^ $: m^^^^^)- 

fcGz;t,|fc|^Af *=i fcez+,|fc|^Af 

We also have to take into consideration the term 

71 

E^A^V^^ne)W = AV(^). 
Finally, we obtain: 

n— 1 n 

-C-^ E ^)Am 

fcGZ+,|fc|s£7V 

The lemma is proved. D 
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